Implementation of universal quantum gates based on nonadiabatic geometric phases 
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We propose an experimentally feasible scheme to achieve quantum computation based on nonadi- 
abatic geometric phase shifts, in which a cyclic geometric phase is used to realize a set of universal 
quantum gates. Physical implementation of this set of gates is designed for Josephson junctions and 
for NMR systems. Interestingly, we find that the nonadiabatic phase shift may be independent of 
the operation time under appropriate controllable conditions. A remarkable feature of the present 
nonadiabatic geometric gates is that there is no intrinsic limitation on the operation time. 
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Quantum computers could efficiently solve certain sig- 
nificant problems which are intractable for classical com- 
puters EJ|. The physical implementation of quantum 
computation (QC) requires a series of accurately control- 
lable quantum gates. These gates may be implemented 
experimentally by using controlled dynamic or geometric 
operations. It is remarkable that geometric operations 
based on adiabatic passages depend only on the global 
feature of the path executed M , and therefore provides 
a possible fault-tolerant way to perform quantum gate 
operations ^-^]. Recently, several schemes for adiabatic 
geometric QC were proposed by using nuclear magnetic 
resonance (NMR) H], trapped ions ||, or superconduct- 
ing nanocircuits |6|]. In particular, an experimental re- 
alization of the conditional adiabatic phase shift was re- 
ported H with the NMR technique. However, the distor- 
tions from the adiabatic approximation were also seen in 
the NMR experiment. Moreover, the adiabatic condition 
requires that the evolution time must be much longer 
than the characteristic time tq of a qubit system; while 
the evolution must be completed within the decoherence 
time. This constraint leads to an intrinsic limitation on 
the operation time of quantum gates and seems to be a se- 
rious obstacle to the physical implementation of some ge- 
ometric QC schemes, especially for those with solid state 
systems where the decoherence time is very short. 

Even for the single qubit gate operations based on the 
adiabatic Berry phase, they might be hard to implement 
with the current experimental technique Therefore, 
it is important to overcome the operation time limitation 
in geometric QC set by the adiabatic condition. A two- 
qubit nonadiabatic geometric phase shift gate was pro- 
posed jio) recently for NMR systems. However, it is still 
an interesting open problem to design a general scheme to 
achieve the universal set of nonadiabatic geometric gates 
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which can be implemented in various physical systems. 

In this Letter, we propose a general scheme to achieve 
the universal set of quantum gates based on cyclic nona- 
diabatic geometric operations. The present scheme is 
experimentally feasible, with the required experimental 
techniques not more stringent than those for dynamical 
gate operations. Similar to the case with adiabatic pas- 
sages, the geometric gates based on nonadiabatic cyclic 
operations also depend only on some global features [|ll| , 
which make them robust to certain computational errors. 
Comparing with a nonadiabatic two-qubit gate proposed 
in Ref.[10], our scheme has at least two distinct advan- 
tages: (i) both a two-qubit gate and two noncommutable 
single qubit gates in a realistic system are designed based 
on a nonadiabatic geometric method, the latter being 
also highly nontrivial and useful in QC || ; (ii) the time- 
independent nonadiabatic phase shift may be realized. 
It is also remarkable that the set of gates designed are 
not only all-geometric but also all-nonadiabatic, and the 
scheme is applicable for several potential physical sys- 
tems. 

For universal QC, we need only to achieve two kinds 
of noncommutable single-qubit gates, and one nontriv- 
ial two-qubit gate |l2). Before the design of geomet- 
ric quantum gates, we show first how to calculate the 
adiabatic and non-adiabatic geometric phases. For a 
spin-1/2 particle subject to an arbitrary magnetic field 
B, the nonadiabatic cyclic Aharonov-Anandan (A A) 
phase is just the solid angle determined by the 
evolution curve in the projective Hilbert space-a unit 
sphere S 2 . Any two-component 'spin' state \ip) — 
[e- lv/2 cos(6/2), e lv/2 sin(9/2)} T may be mapped into a 
unit vector n = (sinOcosip, sinOsintp, cosO) in the pro- 
jective Hilbert space via the relation n = (ip\ a\ ip), 
where T represents the transposition of matrix. By 
changing the magnetic field, the AA phase is given by 
7 = — | f c (l — cos8)dip, where C is along the actual evo- 
lution curve on S 2 , and is determined by the equation: 
<9tn(t) = — /iB(t) x n(t)/h. This 7 phase recovers Berry 
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phase in adiabatic evolution [13]. 

We then show how to achieve the universal set of quan- 
tum gates based on nonadiabatic geometric A A phases. 
The single qubit Hamiltonian H is chosen to go through 
a cyclic evolution with period r in the parameter space 
{B}. We consider the process where a pair of orthogonal 
states \ip±) can evolve cyclically. It is necessary to first 
decide the above cyclic evolution states. A phase differ- 
ence between and can be introduced by cycli- 
cally changing H. The phases acquired in this way would 
contain both a geometric and a dynamical component. 
The dynamical phase accumulated in the whole process 
can be removed by a simple method to be described later 
, and thus only the geometric phase needs to be con- 
sidered at present. By taking into account the cyclic 
condition for \ip±) and removing the dynamical phase, 
we have the relation U(r)\ip±) = exp(±i"/)\ijj±} , where 
U(r) is the evolution operator. Here we have also used 
the result that 7(— n(0)) = —7(11(0)) at any time if the 
two initial states correspond to ±n(0) We now write 

an arbitrary initial state as \ipi) — a_|-|^)+) + a. | "0— } with 

a ± = (V'±I'0)j an d express the two cyclic initial states as 
= cos^\0)+sin^\l) and \ip-) — — szn||0)+cos||l), 
where |0) and |1) constitute the computational basis for 
the qubit. The final state at time r is found to be 
1^/) = U{x,l)\i])i), where 
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It is straightforward to verify that two operations 
J7 1 (Xi,7i) and U 2 (x2,l2) are noncommutable unless 
sinji sin72sin(x2 — Xi) = 0. Since two kinds of noncom- 
mutable operations constitute a universal set of single-bit 
gates, we achieve the universal single-bit gates by choos- 
ing xi 7^ X2 (mod 2tt) for any nontrivial phases 71 and 72. 
For example, the phase-flip gate Ui = exp(— 2i7i|l)(l|) 
(up to an irrelevant overall phase) is achieved at x = 0; 
the gate J7 2 = exp(ij 2 o~x) is obtained at x = 7r /2, which 
produces a spin flip (NOT-operation) when 72 = n/2 
and an equal-weight superposition of spin states when 
72 = 7r/4. 

In terms of the computational basis {|00),|01),|10), 
1 11)}, the unitary operator to describe the two-qubit gate 
is given by U tq = diag(U^o iX o\, C/( 7 i iX i)) under the con- 
dition that the control qubit is far away from the res- 
onance condition for the operation of the target qubit. 
Here 7 5 (x S ) represents the geometric phase (the cyclic 
initial state) of the target qubit as long as the state of the 
control qubit corresponds 8 = 0, 1 (6 denotes the state 
of control qubit). Following Ref. Ea], we find that the 
unitary operator U tq is a nontrivial two-qubit gate if and 
only if 7 1 7^ 7 or x 1 7^ X° (mod 2vr). Therefore, all 
elements of QC may be achievable by using nonadiabatic 
cyclic geometric operations. 



We now describe briefly how to remove the dynamical 
phase |||| . We let \ip + /_) evolve along the time-reversal 
path of the first-period loop during the second period, 
namely, the same loop as before is covered backwards by 
Hb = — jtiB- a /2. This process may be simply realized 
by reversing the effective magnetic field with B(2r — t) = 
— B(i) on the same loop of the first period [0, r). We thus 
have Hb(2t — t) = —H B (t) on the loop. As a result, the 
geometric AA phases accumulated in the two periods will 

add and the dynamical phase will be cancelled. This is 

(2) 

because the dynamical phase 7^ for the second period 
has the same magnitude as that(7^) for the first period 
but with the opposite sign, i.e., 7^ = dt(tp\HB(t)\ip) 

=/ T dt{i/,\H B (2r - *M= - !Z dt{i>\H B (t)\i,)= - 7 «, 
where the fact that | </>(£)) = \ip(2r — t)) and \ip(t)) repre- 
sent the same quantum state is used. 

So far, we have proposed a general scheme to achieve 
nonadiabatic geometric QC. It is important to further 
consider implementing this scheme with real physical sys- 
tems. Here, we illustrate this implementation by two ex- 
amples. The first one is an NMR system j|,[l4]], where 
the Hamiltonian for a single qubit is given by 

H = — — (ujoo- x cosujt + ujQO-ySinujt + wi<7 z ), (2) 

with LUi = nBi/h. Equation (^) can be solved analytically 
]l3[ . In terms of explicit form of the solution n(x,wt) 
represented in Ref. |L3| , it is found that the initial states 
\ip±) with x = arctan[LUo/(uJi + uj)] takes a cyclic evo- 
lution with the period t — 2tt/uj JTo[] , and the evolution 
path is the curve on a unit sphere swept by the unit 
vector ±n(x,ojt). The corresponding geometric phase is 
given by j± = ±7r(l — cosx) for one cycle. Based on the 
noncommutable criterion mentioned before, we may use 
any two processes with different wo/(wi + w) to achieve 
two noncommutable qubit gates. The advantage of the 
above nonadiabatic gates is that there is in principle no 
limitation on the magnitude of to. It needs to be noted 
that, the gates may not be practical by using the 
field B in Eq.(|) since 7 = (tt) as X = (tt/2). This 
problem can be solved by rotating the field. It is seen 
that the parameter x for the initial cyclic state may be 
changed by rotating the symmetric axis for the field (||). 
In the rotated coordinates, B' = R(y, x' ~ x)B (where 
R(y, 9) represents the rotation of angle 8 around the y- 
axis ) and n' (= R(y, x' — x) n (X; w 0) because of the 
spherical symmetry of the system. Thus x may change to 
any required x' f° r implementation of the quantum gate, 
with the geometric phase being unchanged since the area 
swept by n' is the same as that by n. For example, if 
the magnetic field is B' for x' = 7r /2, we may achieve 
the gate U2 — exp(ij20~ x ) with 72 = 27r(l — cosx) (where 
the factor of 2 arises from the evolution of two cycles, 
which is adopted to remove the dynamic phase). A sim- 
ilar method may be employed to achieve the two-qubit 
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operation. The spin-spin interaction in NMR is very well 
approximated by Hj — Ja\a 2 j2. We may use the initial 
state \tjj±) with x S — arctan{uJo / '[uji + (26 — 1) J + ui}} 
to achieve a nonadibatic cyclic two-qubit gate U tq . The 
state of control qubit is (almost) not affected by any op- 
eration of the target qubit if uj\ of the target qubit is 
chosen to be significantly different from u>\ of the control 
qubit (i.e., \oj[ — u>°\ » J). 



Zhu, Wang, Fig .1 



10 20 30 40 50 




1 .2 



50 100 150 200 250 300 

x/x 



FIG. 1. Conditional geometric phases in an NMR system 
for (a) Wj = 0.8J and (b) Wi = J — u with luq = 2\/l5J. 
The solid lines represent actual phases for cyclic evolutions, 
while the dotted ones are calculated under the adiabatic ap- 
proximation. Open circles and squares denote 7 and 7 1 , 
respectively. 

As a typical case, we now show how to achieve 
the conditional geometric phases for the two-qubit gate 
^(2tt 37i72)" ^ S as a f unc ti° n °f T / r o (with To = 2ll I '(jJq) 
is plotted in Fig. la and Fig. Ib, where cuq is set to be 
J. First, we assume that u\ — aJ with < a < 1, 
the phase shifts for Ur$ r 3^/2) are calculated for r = 
2v / 15/[(I - a)J\. The curve 7* ~ r is plotted in Fig. lb 
for a — 0.8(the operation time r = 5Vl5/(2J)). Sec- 
ondly, we choose uj\ — J — u> if u>i is changeable, then the 
phase 7 5 would be independent on the operation time. 
The time-independent conditional phase shifts are clearly 
seen in Fig. lb, namely, the state always evolves along 
the same closed curve in the projected Hilbert space for 



the chosen parameters. This also indicates the geometric 
feature of nonadibatic AA phase. Normally, the state in 
the projected Hilbert space is controlled by varying the 
effective field. 

The second example is provided by the charge qubit 
using Josephson junctions, which has been studied in 
Rcf. H with adiabatic passages. The generalization to 
non-adiabatic quantum computation is important for this 
kind of solid state qubits since the decoherence time in 
these systems is typically short. The single-qubit(Fig.2a) 
consists of a superconducting electron box formed by an 
asymmetric SQUID with the Josephson coupling E\ and 
E2, pierced by a magnetic flux $ and subject to an ap- 
plied gate voltage V x = 2en x /C x with 2en x as the offset 
charge and C x as the capacitance of the electron box. 
In the charging regime (where E\p « E c h with E c h 
as charging energy) and at low temperatures, the system 
behaves as an artificial spin- 1/2 particle in a fictitious 
magnetic field |1| 

B = {Ejcosa, —Ejsina,E c h(l — 2n x )} (3) 



where Ej = y E 2 _ + 4EiE 2 cos 2 (ir&/Q ) 1 tana — 

E-tan(TT<$>/$ )/E + with E± = E\ ± E 2 and $ = h/2e. 
By changing V x and $, Eq.(^|) draws a curve in the pa- 
rameter space {B}. We here study a specific process 
described by 7r$(t) = <& atan[E + tan(tjjt) / EJ\, n e x (t) 
1 — (EjctgxQ + fiu)/E ch ]/2. The fictitious field is a 
rotating field with a constant frequency(Fig.2b), which 
guarantees that the angle xo — arctan[E j / (B z (t) — Tilj)] 
is time-independent. We can find that the state denoted 
by the vector n(xo,— wt) undergoes a cyclic evolution 
with the period r = 2it/u Therefore, we can obtain 
a unitary operator ([l]) in the charge-qubit system, where 
X = Xo an d the AA phase 7 is the half solid angle swept 
by the vector n(xo, —cot). 

The interaction between different charge qubits may 
be realized by coupling two asymmetric SQUIDS ca- 
pacitively. If the coupling capacitance Ci is smaller 
than the others, the field on the target qubit is given 
by Eq.(||), but the z-component is replaced by B & z = 
E c h(l — 2n%) + Ei(n x c — S), where n% c represents the 
offset charge in the control qubit and Ej is the coupling 
energy || . Obviously, the j s phase of the target qubit in 
the decoupled case is different from 7 even though ($, n%) 
varies in the same way. If the offset charge n x c is time- 
independent in the process , the state n(x 5 ,— tot) with 
X s = arctan[Ej / (B 5 Z — hui)) still undergoes a cyclic evo- 
lution. Thus the two-qubit operator U tq may be obtained 
similar to the case in NMR. The elimination of the adia- 
batic condition for two-bit geometric gates is significant 
since the coherence time for two qubit gates is typically 
much shorter than that for single qubit gates. 
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in NMR systems addressed before. 

Finally, we emphasize that the advantages of our 
scheme are: i.) the designed quantum gates are universal 
and can handle arbitrary QC without the intrinsic limita- 
tion on operation time; ii.) the scheme is essentially all- 
geometric and is robust to certain computational errors; 
iii.) the physical implementation of the scheme can now 
be realized by realistic Josephson junctions and NMR 
systems. Nevertheless, it seems to be a limitation of the 
method that the nonadiabatic phase is more sensitive to 
fluctuation in the trajectory area than that of the adia- 
batic phase. 
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FIG. 2. (a) A schematic Josephson charge qubit(Ref. [6]). 
(b) The fictitious magnetic field (in unit of fj,ev) in 
the designed cyclic process for nonadiabatic geometric 
gates. (c) Geometric phases to achieve a NOT gate 
for E 2 = 4Ej = 6.25fiev, E ch = 5.0(E L + E 2 ), and 
Xo = arccos{3/4). The inset shows the phase shift for the 
Hadamard operation as xo = arccos(7/8). The solid lines 
represent the actual phases(7) for cyclic evolutions, while the 
dotted ones(7 a ) are calculated under the adiabatic approxi- 
mation. 

To see clearly the advantages of nonadiabatic geomet- 
ric qubit gates, we compare the operation time between 
nonadiabatic gates and adiabatic gates, with the phase 
shift 7r/2 (7r/4) corresponding to a NOT (Hadamard) op- 
eration. The geometric phases for single-bit operations 
are shown in Fig. 2c with certain parameters. It is seen 
from the main panel of Fig. 2c that only if r is longer 
than 70t (= Ti/Ej) , the adiabatic phase j a deviates 
from the actual nonadiabatic phase 7 within 10% error. 
The similar results can be obtained for other parame- 
ters^. g., see the inset of Fig. 2c for 7 = 7r/4). Note 
that the coherence time achieved by current technology 
is 30 ~ 40to H and thus the theoretical proposal based 
on the adiabatic phase is not accurate enough to achieve 
this kind of quantum gate experimentally. In addition, 
Fig. 2c also clearly shows that nonadiabatic phases can 
be independent of the operation time, similar to the case 
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